clear;
T=2;
N=1000;
dt=T/N;%diataT
t=0:dt:T;%时间
y0=1;%初始条件
y1=zeros(1,length(t));%建立一个全是0的向量
y1(1)=y0;%将初始条件放入
for n=1:length(t)-1
    %y(n+1)=y(n)+y(n)*dt;
    y(n+1)=y(n)+(sin(y(n))*t(n)+1)*dt;
end
plot(t,y);
%利用holdon更改划分，判断怎么样的精度合适
hold on